Optimal degeneracy for the signal-recycling cavity in advanced LIGO 
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As currently designed, the signal-recycling cavity (SRC) in the Advanced-LIGO interferometer is 
degenerate. In such a degenerate cavity, the phase fronts of optical fields become badly distorted 
when the mirror shapes are slightly deformed due to mirror figure error and/or thermal aberration, 
and this causes significant loss of the signal power and the signal-to-noise ratio (SNR) of a gravita- 
tional wave event. Through a numerical modal simulation of the optical fields in a simplified model 
of an Advanced-LIGO interferometer, We investigate the loss of the SNR and the behavior of both 
the carrier and signal optical fields, with the SRC at various levels of degeneracy. We show that the 
loss of the SNR is severe with a degenerate SRC, and a nondegenerate SRC can be used to solve 
this problem. We identify the optimal level of degeneracy for the SRC, which is achieved when the 
cavity Gouy phase is between 0.2 and 1.3 radians. We also discuss possible alternative designs of 
the SRC to achieve this optimal degeneracy. 

PACS numbers: 



I. INTRODUCTION AND SUMMARY 

Advanced LIGO Q entails, among other upgrades 
from initial LIGO Q , introducing a signal-recycling mir- 
ror (SRM) at the dark port output of the interferometer 
(see Fig.JTJ. 

The SRM forms the signal-recycling cavity (SRC) with 
the input test mass (ITM), and the SRC and the arm 
cavity (AC) form a coupled resonant cavity. The reso- 
nant property of this coupled cavity can be controlled 
by two parameters of the SRM (position and reflectivity) 
Pi . With different choices of these parameters, the in- 
terferometer can operate in either a broadband, resonant- 
sideband- extraction (RSE) configuration or a nar- 
rowband configuration. The Advanced-LIGO baseline de- 
sign adopts the RSE broadband configuration, with the 
possibility, later, of changing the SRM parameters so as 
to alter the detector noise spectrum, optimizingits de- 
tection of GWs with specific frequency features [5j . 

Signal recycling can also circumvent the standard 
quantum limit (SQLjfor free test masses by altering the 
test- mass dynamics 0. 

However, there is a potential problem in the current 
design of the SRC. 

The SRC and the power-recycling cavity (PRC) are 
both nearly degenerate, and it is well known that a degen- 
erate cavity is not selective for transverse optical modes. 
As a result, perturbations of the cavity geometry will 
cause strong mode coupling Specifically, in initial- 
LIGO and Advanced-LIGO interferometers, figure error 
and thermal aberration of the mirrors (PRM, SRM and 
ITMs) will cause strong optical mode coupling, which 
transfers light power from the fundamental TEMOO mode 
to higher-order modes (HOMs) and reduces the ampli- 
tudes of the radio frequency (RF) sidebands (in both the 
PRC and SRC) and the signal sideband (in the SRC). 

This consequence of the high PRC degeneracy is well 
known: strong mode mixing of the RF sideband has been 
observed in initial LIGO. This problem was so severe that 
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FIG. 1: Illustration of an Advanced-LIGO interferometer. A 
signal-recycling mirror is placed at the dark output port of 
the LIGO interferometer, forming an SRC with the ITMs. 
In this plot we introduced symbols and abbreviations for the 
mirrors and the electric fields at various positions that will be 
used in this chapter. The two ACs are named "online" and 
"offline" for convenience. 



the interferometer had to be operated with lower circu- 
lating power to reduce thermal aberration of the mirrors. 
Measures had been taken to fix the problem, including 
introducing a thermal compensation system (TCS) to ac- 
tively correct the deformations of the mirrors fTol | and re- 
placing bad optical elements that had unexpectedly high 
absorption. However, there is a worry about Advanced 
LIGO, where much higher optical power in the AC will 
cause even worse thermal aberration on the ITMs that 
the TCS might not be able to correct. Miillcr and Wise 
have suggested reducing the PRC degeneracy by mov- 
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ing the mode-matching telescope (MMT) into the recy- 
cling cavities to reduce the beam waist size 0| , and they 
are currently working on practical issues in implementing 
this modified topology in Advanced LIGO [I3 ■ 

The consequence of SRC degeneracy, by contrast with 
PRC degeneracy, have not been clearly investigated. 
Since the GW signal sideband light entering the SRC 
has different resonance conditions from the carrier light 
or the RF sideband light, the three must be investigated 
individually. Specifically, the signal sideband is resonant 
in the coupled two-cavity system formed by the AC and 
the SRC, while the carrier and the RF sideband, roughly 
speaking, are resonant only in the nondegenerate AC and 
the degenerate PRC respectively. It seems that mode 
mixing in the signal sideband is at a level somewhere 
between those of the carrier and the RF sideband. 

In Section IV J of Ref. Thorne estimated the 

strength of mode mixing in the signal sideband using 
the approximation that light propagation in the degen- 
erate SRC is well described by geometric optics. (Be- 
cause a very degenerate cavity accommodates most the 
optical cigenmodes up to very high orders, "light rays" 
with sharp edges are eigenmodes of the cavity as well). 
Thorne found (Eq. (4.54) of Ref. 43]) that for the signal- 
to-noise ratio (SNR) to be reduce by less than 1% due 
to mode mixing , the peak-to-valley mirror figure er- 
ror in the central region (region enclosing 95% of the 
light power) in the SRC has to be less than 2 nm for 
the broadband Advanced-LIGO baseline design, and less 
than lnm for a narrowband configuration. This is in- 
dependent of whether a Gaussian beam or Mesa beam 
[l3j is used, since a degenerate cavity does not distin- 
guish optical modes. These are severe constraints that 
are difficult to achieve with current technology. 

In this chapter, we investigate this SRC degeneracy 
problem using a modc-decomposition-based numerical 
simulation of light propagation, including both the car- 
rier and the signal sidebands, in a simplified model of 
an Advanced-LIGO interferometer. In our simulations, 
we focus on the consequences of phase-front distortions 
of the light in recycling cavities with various levels of 
degeneracy. We make idealizations and approximations 
to simplify the analysis of the interferometer, so long as 
they do not make errors larger than a factor of order 
two. All assumptions and approximations are discussed 
in detail in this chapter and we give a full list of them 
in Appendix iBl Among the most important ones are the 
following: 

(i) We consider only the five lowest-order HOMs and 
focus mostly on the two lowest-order modes, i.e., the 
modes excited by errors in mirror curvature radii. 

(ii) We assume that the sizes of mirrors in the interfer- 
ometers are all large enough for diffraction losses to be 
negligible. 



Assuming that shot noise is the dominant noise source 



(hi) We assume that the interferometer noise is domi- 
nated by photon shot noise, and we ignore radiation pres- 
sure noise. 

(iv) We assume all degrees of freedom of the interfer- 
ometer (cavity tuning, optical element alignment, etc.) 
are fixed to their ideal values as if there were no mirror 
deformation, except for the tuning of the PRC and the 
AC. The PRC and AC are tuned to maximize the total 
carrier light power (sum of power in all optical modes) 
in the AC, for fixed input light power. 

We use our simulations to study the loss of the SNR 
due to phase-front distortions caused by mirror deforma- 
tions, and reach the following conclusions: 

(i) With the degenerate SRC in the current Advanced- 
LIGO baseline design, if we require the loss of the SNR 
due to mirror deformations to be smaller than 1%, then 
the constraint on mirror deformations is severe. In the 
broadband Advanced-LIGO design, if the only type of 
mirror deformation is curvature radius error, then this 
error must be smaller than 2.5m~7m on the ITMs, which 
corresponds to a peak-to-vally figure error smaller than 
lnm~3nm in the central region of the mirrors. In real- 
ity, the mirror deformation is formed by the combination 
of many spatial modes, and when we consider the next 
lowest order spatial modes, as depicted in Eq. I|27|) and 
Fig. [2 we get a constraint ~ 4nm on the SRM. In the 
narrowband design, the constraint on the ITMs is tight- 
ened to ~ 0.4nm. These results are consistent with the 
estimates in Ref. [l3| . 

(ii) When we change the degeneracy of the SRC while 
keeping tunings of other cavities the same, we find for 
Gaussian beams that the loss of the SNR due to mirror 
deformations is minimized when the one-way Guoy phase 
inside the SRC is in the range 0.2^1.3 radians and is 
chosen to be away from HOM resonant peaks (sec Sec lIVjl . 

(iii) We find that it is not practical to add a single 
lens into the SRC so as to focus the beam and reduce 
the cavity degeneracy to its optimal value, because the 
beam must be focused so strongly that the beam size on 
the SRM is of order 10~ 5 m, and the power density on it 
exceeds 10GW/m 2 (This fact was pointed out long ago 
by Bochner Q). 

(iv) We discuss two alternative designs of the SRC for 
achieving the optimal degeneracy. The first is the MMT 
design Muller and Wise suggested originally for reducing 
the degeneracy of the PRC [Uj. In the MMT design, the 
beam size is brought down to ~ 10 m by two mirrors 
so as to achieve the optimal degeneracy. The second is a 
long SRC design, in which the beam size is kept at ~6cm, 
while the SRC is extended to ~4 kilometers to achieve 
the optimal degeneracy. The 4km long SRC design has 
been suggested for various reasons 0, fig , and in practice 
it is possible to fit the SRC into the existing LIGO AC 
beam tubes. 

This chapter is organized as follows. In Sec. [HJ we 
give a brief overview of the mode decomposition formal- 
ism and Hermite-Gaussian modes |l6j . and interpret the 
cavity degeneracy from a modal-space point of view. In 
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Sec. IIIII we describe the Advanced-LIGO interferometer 
model that is used in our simulations. In Sec. II VI we 
summarize the numerical results that come out of our 
simulations, including the constraints on mirror figure 
error and thermal aberration, and the optimal SRC de- 
generacy. In Sec. we discuss the various designs for 
achieving the optimal SRC degeneracy. In Sec. IVII we 
summarize our conclusions. 



II. MODE DECOMPOSITION FORMALISM 
A. Modal decomposition in general 

The mode decomposition formalism for calculating op- 
tical fields in a perturbed interferometer is discussed in 
detail by Hefetz et al in Sec. 2 of Ref.01- We will review 
the general idea briefly in this section. 

One can generally expand the electromagnetic (EM) 
field of a light beam as a superposition of orthonormal 
optical modes: 



E(x,y,z) 



a n U n (x,y,z) . 



(1) 



Though the basis modes U n (x,y,z) are arbitrary in 
principle, it is preferable to use the eigenmodes of the 
cavities of an ideal interferometer, e.g., (i) Hcrmite- 
Gaussian modes, which are eigenmodes of the cavities 
in the Advanced-LIGO baseline design, formed by spher- 
ical mirrors (assuming infinite mirror size); or (ii) Mesa- 
beam modes 01 > which are eigenmodes suggested for 
Advanced LIGO to reduce thermal noise. The complex 
vector space formed by U n (x,y, z) is called the modal 
space, and the EM field in modal space is represented by 
a complex vector a n . (In Sec. IIIII and Appendix IO we 
keep using E to denote this vector.) In the modal space, 
the optical fields of a perturbed interferometer can be cal- 
culated from the unperturbed fields using linear algebra 
only, without numerically solving the wave equation. 

The propagation of the optical field can be described 
by matrices in this modal space. In Cartesian coordinates 
where the z-axis is along the optical axis and the x- and 
y-axis arc transverse, an operator M(x, y, z 2 , z\) trans- 
forms the EM field at position z\ to the field at position 

E(x 2 ,y 2 , z 2 ) = M(x2,y2;xi,yr,Z2, z\) ® E{x\,yi,zi) . 

(2) 

The representation of M in the modal space is given by 



M mn (z 2 , Z ^ = JJJJ U m( x 2,V2, Z 2 ) 

x M(x 2l y 2 ;xi,yi; z 2l zi)U n (xi,yi, zi)dxidyidx 2 dy 2 . 

(3) 

It is convenient to separate these operators into prop- 
agation operators in free space and interaction operators 



describing how the EM fields transform when interacting 
with optical elements. The free-space propagator is given 
by: 



-ik(z 2 —z 1 ) girjn 



(4) 



where k is the wave number and r\ n is the diffrac- 
tion phase associated with the nth optical mode (i.e., 
the extra phase accumulated during propagation besides 
k(z 2 — zi), due to diffraction effects, e.g., the Guoy phase 
of Gaussian beam). 

To write out the interaction operator for an optical 
element, we choose, near the element's surfaces, reference 
surfaces that match the phase fronts of the unperturbed 
eigenmodes. The operator can then be written in the 
general form 



Mr, 



(m\M(x,y)\n) 



(m\e- ikZ( - x ' y) \n) 



(5) 



where Z(x,y) is the optical path light travels from the 
element's entrance reference surface to the element's exit 
reference surface. For ideal optical elements in our model, 
Z(x, y) = constant, i.e., the mirrors or lens exactly match 
the optical modes and there is no coupling between opti- 
cal modes when the light beam interacts with the optical 
elements. For perturbed elements, e.g., slightly deformed 
mirrors due to figure error and/or thermal aberration, 
Z(x,y) is not constant and the optical modes couple to 
each other. Z[x,y) contains contributions from both the 
figure error and the change of refraction index in the ma- 
terial, and is referred to as the distortion function] it 
can be complex, when used to describe lossy optical el- 
ements. This interaction operator must be accompanied 
by the optical element's scalar reflectivity and transmis- 
sivity coefficients to give the true transform of the fields. 

In writing the interaction operators, we adopted the 
short-distance approximation |l8j , where propagation in- 
side optical elements between reference planes is approx- 
imated by a simple non-uniform phase factor kZ{x,y). 
The spatially variable phase error caused by this approx- 
imation, in addition to a factor of the order unity de- 
termined by the geometry of the unperturbed cavity, is 
derived in Sec. 2 of Ref. [Tg| ; its magnitude is 



(0) 



where A is the wavelength of the light and L is the length 
of the cavity. In initial LIGO or Advanced LIGO, L is 
at least ~ 10m, and the phase distortion error is thus 
smaller than 10~ 8 . Since a lnm mirror figure error gives 
a phase distortion of 6 x 10 -4 , A$ is always negligible. 



B. Hermite-Gaussian modes 

In this section, we review briefly the Hermite-Gaussian 
modes (see Chapters 16 and 17 of Ref. ^||) that are used 
as basis modes in our simulation. 
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U n (x,z) = - 



A Hcrmitc-Gaussian mode of beam waist size wq is 
given by 

2\ V V 1 V 7 * fV2x\ 



(7) 



7r J \2 n n\w(z) J \ w(z) 

1 ik 



xexp — x 



w{z) 2 2R(z), 



xexp [ i ( n+ - ) r?(z) ) , 



where H n (x) is the Hermite polynomial and R(z), w(z), 
and ?;(z) are the curvature radius of the phase front, the 
beam spot size and the Guoy phase, respectively, given 
in terms of the Rayleigh length zq = ttwq/X by 



R(z) 



z + 



and 



w(z) 



r]{z) 




tan 



(8) 



These Hermite-Gaussian modes are exact solutions to the 
paraxial wave equation in one dimension, and they form a 
complete orthogonal basis in the solution space. We use 
Hermite-Gaussian modes as basis modes in the modal 
space 



E(x,y,z) 



U m (x, z)U n (y, z)exp(-ikz) , (9) 



where each transverse mode is labeled by two integers 
(m, n) corresponding to the directions x and y. In this 
chapter, we consider only the five lowest-order symmetric 
HOMs, i.e., modes with even m and n, and m+n < 4. We 
omit the two m + n = 1 modes because they can be cor- 
rected by the tilt control system. We omit modes with 
m + n > 4 to reduce computational cost, and because 
modes with higher orders are suppressed by stronger 
diffraction loss, their coupling to the fundamental mode 
is weaker when mirror deformations are smooth. More- 
over, mode-coupling behaviors due to mirror deforma- 
tions are qualitatively the same for all modes, and there 
is no need to include more modes with order m + n > 4 
for an order of magnitude estimation. We also focus our 
attention on the second order (m + n = 2) modes since 
they are most likely to be the dominant deformations 
present in Advanced LIGO due to mirror figure error or 
thermal lensing, and thus are the lowest order perturba- 
tions of Hcrmitc-Gaussian modes after mirror tilts have 
been suppressed by control systems. 

The propagation operators for Hcrmitc-Gaussian 
modes arc given by Eq. Q, with the diffraction phases 
replaced by the Gouy phases of the (m, n) Hcrmitc- 
Gaussian modes: (m + n + l)r)(z). The interaction op- 
erators defined in Eq. (JSJ are derived analytically for 
Hermite-Gaussian modes in Section 2 of Ref . , assum- 
ing that the mirror radius is much larger than the beam 
size. For the mirror deformations considered in this chap- 
ter, this approximation produces a few percent errors in 



the coupling strengths between optical modes (measured 
by components of the interaction operators). In fact, 
these errors are less than 10 -3 for the coupling between 
the fundamental mode and the five HOMs we considered, 
which is the leading order effect we would like to investi- 
gate in this chapter. Only for coupling between HOMs, 
which is a higher-order effect in changing the SNR, does 
the error caused by the finite size of mirrors become as 
large as 1 ~ 10%. Also, because Gaussian modes are 
cigenmodes of cavities formed by spherical mirrors with 
infinite size, it is in fact self-consistent that we treat the 
mirrors as large when using Gaussian modes as eigen- 
modes. As long as we use eigenmodes consistent with 
mirror size, even if we use infinite-size mirrors, our esti- 
mates are valid up to fractional errors of order the differ- 
ence between the Hermite-Gaussian modes and the true 
cigenmodes of the cavity with finite-size mirrors. Thus, 
we approximate the mirrors as having infinite sizes and 
use Gaussian modes throughout the chapter. 

Gaussian beams have spherical phase fronts, and are 
thus eigenmodes of cavities formed by spherical mir- 
rors. There are simple relations between the geometry 
of the cavity (measured by, e.g., cavity g-factor, and 
mirror radii of curvature) and the physical properties of 
the Gaussian eigenmodes (e.g., Guoy phase, beam waist 
size, waist position) that are available in Chapter 19 of 
Ref. [3. 

A very useful formula relating a cavity's one-way Guoy 
phase r) (the Guoy phase of the cavity's fundamental 
cigenmode) to the cavity g-factor is: 



77 = arccos ^fg . 



(10) 



The cavity g-factor is the product of the two mirror g- 
factors. If the curvature radii of the two mirrors are R\ 
and i?2, and the length of the cavity is L, then the mirror 
and cavity g-factors are defined as 

9i = 1 - -5- = 1, 2) and g = g x g 2 ■ (11) 
Hi 

The cavity is stable (i.e., the cavity geometry supports 
Hermite-Gaussian modes as eigenmodes) if and only if 
0<9 <1. 

From Eq. QJ, we see that it is the Guoy phase that 
distinguishes optical modes (modes with different m + 
n for Hcrmitc-Gaussian modes). If the one-way cavity 
Gouy phases are very close to or w, the round trip phase 
shifts in the cavity are almost the same for all modes, so 
if one of the optical modes is tuned to be resonant, so are 
the others. Such a cavity is thus called degenerate. In 
terms of g-factor, a cavity is degenerate if g is very close 
to or 1. 

To conclude this section, we use the baseline design of 
the AC and the PRC in Advanced LIGO as examples to 
demonstrate quantitatively the degeneracy level of cavi- 
ties. The baseline curvature radii of the PR mirror and 
the test masses are 0] 



-Retm = Rxtm = 2076.4m, 
R PR = 1194.7m, 



-R1TM2 



-1186.4m, 
(12) 
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where Rjtm and i?iTM2 are curvature radii of the ITM 
seen from inside the AC and the PRC, respectively, and 
the ITM is convex as seen from the PRC. The cavity 
lengths are d A c = 4000m and g?prc = 8.34m, and the 
Rayleigh lengths and Guoy phases are 

zoac = 390.9m, VAC = 0.39, (13) 
^oprc = 82.1m, 77PRC = 4.9 x HT 4 . (14) 

In the AC, the Rayleigh length is clearly much shorter 
than the typical distance 4km that carrier light travels 
in the cavity, i.e., the light propagation is in the strong- 
diffraction zone, which indicates a nondegcncrate cavity. 
More rigorously, the Guoy phase, corresponding to a fre- 
quency shift of 

Az/ = — ^— 77AC = 4.6kHz , (15) 

is much larger than the bandwidth of the AC (~ 15Hz). 
This means, the Guoy phase breaks the degeneracy be- 
tween the Gaussian modes with different orders (differ- 
ent m + n), i.e., when the cavity is tuned to have the 
fundamental mode in resonance, nearly all other HOMs 
are suppressed. Of course there are always HOMs with 
round trip Guoy phases (mode 2ir) close to that of the 
fundamental mode by coincidence, but these HOMs that 
resonate with the fundamental mode generally have very 
high orders (except for very bad choices of the Guoy 
phase) and are thus strongly suppressed by diffraction 
losses. 

In the PRC, the Rayleigh length is longer than the 
length of the cavity, but shorter than the typical distance 
RF sideband light travels inside the cavity after we take 
count of the number of round trips (~ 50). Therefore, the 
RF sideband propagation in the PRC is still in its strong 
diffraction zone. However, the frequency shift of Au = 
2.8kHz is much smaller than the bandwidth of the PRC 
(~ 100kHz), and consequently the Guoy phases are close 
to and it (mode 2tt). The PRC, therefore, although 
not degenerate to the extreme level that geometric optics 
becomes valid, accommodates tens of low-order HOMs 
together with the fundamental mode, and is thus highly 
degenerate. 

Finally, through the example above, we can see that for 
Hermite- Gaussian modes, reducing the beam waist size 
wq (i.e., reducing Rayleigh length z$) and/or increasing 
the cavity length will reduce the degeneracy of the cavity. 

III. ADVANCED-LIGO INTERFEROMETER 
MODELING 

In this section, we describe our simplified model of an 
Advanced-LIGO interferometer, and the way our simula- 
tions work. 

In our simulations, we study an Advanced-LIGO inter- 
ferometer in equilibrium with static optical fields. We use 
the standard Advanced-LIGO optical topology displayed 



in Fig. 2] and the input light is a pure (0,0) Hcrmitc- 
Gaussian mode coming in from the PR mirror. We con- 
sider both broadband and narrowband interferometer de- 
signs. The interferometer parameters for the broadband 
detector are chosen as their values for the Advanced- 
LIGO baseline design. The par ameters in the two de- 
signs are listed below [Til HH, Il9l| , where we begin using 
the following subscripts to denote different mirrors and 
cavities throughout the whole chapter (see Appendix [S] 
for a list of symbols and subscripts): "bs," "i," "e," "p," 
and "s" stand for the beam splitter (BS), ITM, ETM, 
PRM, and SRM; "ac," "pre," and "sre" stand for AC, 
PRC, and SRC. 

(i) Cavity macroscopic length: The ACs both have 
L = 4000m; the lengths between the PRM and the two 
ITMs are denoted l\ and 1% and referred to as Michclson 
lengths; the lengths between the SRM and the two ITMs 
are denoted I3 and I4. It is convenient to define common 
and differential lengths: 

Zprc = - ^ - = 8.34m , Z src = 3 ^ 4 = 8.327m , 

(16) 

and 



In our model, there is no macroscopic asymmetry be- 
tween the two Michclson arms. Therefore only a micro- 
scopic tuning value is assigned to d in the next paragraph. 

(ii) Cavity microscopic tuning: The carrier light gets 
the following phase shifts during a single trip in the AC, 
PRC, and SRC: 

ac = P c = , 4>l c = 0.06 , Cc = * - 1-556 , 

(18) 

where the superscripts on the SRC phasing denote broad- 
band (B) and narrowband (N). An asymmetry in the 
Micheleson arm lengths is introduced because, among 
other reasons, we choose the homodyne readout scheme 
where a tiny amount of the carrier light power goes to- 
ward the dark port and beats with the resonant signal 
sideband to give the detector output. For this purpose 
solely, the asymmetry is specified at the microscopic level 
as the phase difference the carrier accumulates in the two 
Micheleson arms: A(f> = oj^d/c = 0.01, so that about 1W 
of carrier power goes into the SRC. 

(iii) Mirror power transmissivity: 

tf= 0.5%, t%= 76ppm, ^ = 5.9%, 
tf = 7% , if 2 = 0.3% . (19) 

We assume lossless mirrors throughout our simulation, 
so the amplitude reflectivity and transmissivity are com- 
pletely determined. 

(iv) Mirror curvature radii: These are given in Eq. 1|12|) 
except for the SRM. In this chapter, we will change the 
SRC degeneracy through the value of its Guoy phase, and 
assume that the geometry of the SRM always matches 
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our choice of degeneracy. The corresponding SRM cur- 
vature radius is relevant to nothing but the SRC Guoy 
phase, so we do not specify it explicitly 

The differences between the broadband and the nar- 
rowband designs are all in the choice of the SRM trans- 
missivity and the SRC tuning, as we mentioned in Sec. [I] 
The complex optical-sideband resonant frequency in the 
coupled SRC and AC two-cavity system is given by Eq. 
(13) of Ref. [U: 

ic n + r p e 2 ^=" . 

where X g and e arc positive and are the resonant fre- 
quency and the decay time. With our choice of signal- 
recycling parameters above, we get resonant sideband fre- 
quencies X 9 /2tt = 228Hz and \ g /2ir = 1005Hz. The 
actual resonant frequency is luq — A g , i.e., the down- 
converted signal sideband. Depending on the sign of the 
SRC detuning, the interferometer response is of interest 
for only one of the two sidebands. 

The PRM transmissivity and the PRC tuning are cho- 
sen such that the PRM is impedance matched to the 
AC and the total carrier power (summed over all opti- 
cal modes) in the AC is maximized for fixed input light 
power. 

The only interferometer control we do in simulations 
when mirrors arc slightly deformed is to optimize the to- 
tal carrier light power in the AC by adjusting the tuning 
of the PRC and the AC. To optimize the carrier light 
power in the AC, instead of modeling the control signal, 
we look directly at the power at each equilibrium state, 
i.e., we do a static pseudo-control and do not model the 
dynamical response of the interferometer during the con- 
trol process. With our choice of interferometer parame- 
ters given above, and a 125W input light power, the car- 
rier light power in an ideal interferometer is ~ 825kW, 
the power recycling factor is ~ 18, and the carrier light 
going toward the dark port is ~ 1W. 

Since we are interested in the reduction of the signal 
sideband power due to mirror deformations, we consider 
in our simulations only the carrier light and the down- 
converted signal sideband, and omit the RF sideband and 
other sidebands used for control purposes, as well as the 
up-convcrtcd signal sideband, since the interferometer is 
tuned to be sensitive only to the down-converted signal 
sideband with the best sensitivity (assuming white noise) 
at frequency /o — A g /2-7r. The distortion of the phase 
fronts also affect the control sidebands, but we are not 
concerned with several percent loss of the SNR of the con- 
trol signals. However, there is a problem associated with 
tilt control signals, whose SNRs are proportional to the 
amplitude of the (1,0) and (0,1) Hermite-Gaussian modes 
excited by mirror tilts. If we choose to suppress HOMs 
for the signal sideband by using nondegenerate recycling 
cavities, the SNR of the tilt control signals entering the 
recycling cavities will also be strongly reduced. Corre- 
spondingly, choosing the degeneracy level of the SRC 
will entail compromises between inputs on the signal and 



control sidebands. This problem is not considered in this 
chapter, and is left for future investigation. 

After setting up our model of an ideal interferome- 
ter, we introduce perturbations to its mirrors. We do 
not model thermal lensing of the mirrors. We assume 
that all tilts and misalignments of the mirrors are cor- 
rected by the control system. Because of the high com- 
putational cost associated with our mode decomposition 
method, we limit our figure errors to simple profiles so 
that they generate, at the leading order, coupling be- 
tween only a few optical modes (< 6). One focus of our 
study is mirror curvature radius errors, which is the most 
interesting type of deformation, since it can be generated 
effectively by thermal lensing of the ITMs. We assume 
that the beamsplitter is perfect, because deformations of 
the beamsplitter introduce complicated phase-front dis- 
tortions which have no qualitative difference from those 
introduced by other mirrors. 

At the interferometer output, we assume that there is a 
mode cleaner that filters out all HOMs in the carrier and 
signal sidebands, so that the shot noise is proportional 
to the square root of the output carrier power in the 
fundamental mode [i.e., (0,0) Gaussian mode]: 




Here the superscript "C" stands for carrier, and the sub- 
script labels the mode. The signal power comes from 
beating the signal sideband against the carrier, both tak- 
ing only the fundamental mode, so 

where the superscript "S" stands for signal sideband. As- 
suming shot noise dominates, we have 

SNRcx-y/i^, (23) 

i.e., the SNR is directly proportional to the signal side- 
band amplitude in the fundamental mode. 

When we take into account radiation-pressure noise, 
the change of the SNR becomes more delicate. However, 
since the radiation-pressure noise is determined by the 
total light power on the test masses, it is presumably 
less sensitive to the mode structure of the light. When 
radiation- pressure noise is important, although the loss 
of the SNR is not given by Eq. 12/il) , that equation is still 
a rough measure of the loss of the SNR. This argument 
applies to all other non-optical noise sources. Radiation- 
pressure noise is thus omitted in our simplified model, 
and will be studied in more sophisticated future simula- 
tions pn . 

In the interferometer model described above, we cal- 
culate the signal sideband in two steps. In the first step, 
we propagate the input carrier light (Nd:YAG laser) with 
frequency fo = 2.82 x 10 14 Hz through the interferome- 
ter to build up the static carrier-light field. In the sec- 
ond step, we assume a sinusoidal gravitational wave of 
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frequency f g propagating perpendicular to the detector 
plane with only "+" polarization, i.e., effectively, it dif- 
ferentially shakes the ETMs sinusoidally with frequency 
f g . To leading order in the GW strain, two signal side- 
bands of frequencies fo ± f g are generated at the ETMs 
with exactly the same mode structures as the carrier 
field there. We propagate the down-converted sideband 
through the interferometer to build up the static signal 
sideband field. Repeating this second step with various 
GW frequencies, we map out the frequency response of 
the detector. Repeating both steps with various SRC 
geometries and levels of degeneracy, we can study the ef- 
fect of the recycling-cavity degeneracy on the influence 
of mirror deformations on signal response and the inter- 
ferometer's noise spectrum. 

The EM field in the interferometer (a system of coupled 
optical cavities) can be written formally as 



E = E, 



pump 



Pr.t.E. 



(24) 



Here -E pU mp is the pumping field that contributes di- 
rectly to the E field, and P r .t. is the round-trip prop- 
agator, which consists of free propagation operators and 
interaction operators that describe the propagation of E 
through the interferometer and back to itself. The spe- 
cific forms of -Epump aud P r .t. depend on location in the 
interferometer. For example, the field E c i T i in Fig.Q] i.e., 
the circulating field in the online AC at the ITM going 
toward the ETM, can be written as 



-E'en! — tiTiiE in i + r i r c MiiP &c iM c iP as .iE C ni , 



(25) 



where, e.g., Mn and T;i are reflection and transmission 
operators of the ITM in the online cavity, and we use 
subscripts "1" and "2" to denote the online and offline 
cavities. We can write out a set of coupled equations in 
the form of Eq. (|2*4*jl for all fields labeled in Fig. ^ and 
solve them numerically by iteration, as has been done in 
the FFT simulation code for optical fields in LIGO [l4j . 
In principle, we can also solve for each field in terms of 
the input field E- ln o by directly taking the inverse of all 
operators of the form (I— P r .t.) {I is the identity matirx). 
In the FFT code, hundreds of modes are included, and it 
is computationally difficult to take the inverse of all the 
large matrices, which are sometimes nearly singular. In 
our simulation however, as we consider only six modes, 
we find it more efficient to directly invert the matrices 
instead of iterating the fields. 

In Appendix [CJ we write out the field coupling equa- 
tions explicitly and solve for all the carrier and signal 
sideband fields. 



IV. MIRROR FIGURE ERROR AND OPTIMAL 
DEGENERACY 

In this section, based on results of our simulations of 
the simplified Advanced-LIGO model set up in Sec. II I II 
we try to answer the question how much loss of the SNR 



is caused by mirror deformations with various spatial 
modes and magnitudes, and how does this loss of the SNR 
depend on the degeneracy of the SR cavity? 

For mirror deformations, we consider mostly mirror 
curvature-radius errors. At leading order, we need only 
consider the (2,0) and (0,2) Hermite-Gaussian modes ex- 
cited by curvature radius error. 

We also consider a case in which, besides the curva- 
ture radius error, there is deformation of the SRM with 
higher-order spatial modes. We are interested in this 
case since, if the SRC is nondegenerate, the optical eigen- 
modes of the SRC and the coupled SRC- AC cavities un- 
der the mirror deformation are different from those of the 
coupled PRC- AC cavity, and there might be a substantial 
loss of the SNR due to the mode mismatch between the 
carrier light (mostly in the PRC- AC cavity) and the sig- 
nal sideband light (mostly in the SRC- AC cavity). Thus 
we would especially like to see how the loss of the SNR 
depends on the degeneracy of the SR cavity in this case. 

For signal sidebands, we consider mostly the reso- 
nant signal sideband with frequency fo — ^g/2iT given 
by Eq. H2U|) . and only at the end of this section do we 
consider signal sidebands with frequencies varying from 
50Hz to 1000Hz. 



A. Curvature radius error on the ITMs: 
Broadband configuration 

First, we consider a broadband interferometer with 
curvature radius errors on the ITMs that simulate the 
thermal lensing effect, and we assume all other mirrors 
are perfect. The curvature radii of the ITMs (i?iTM = 
2076.4m) are changed by Ai?iTM = 5m, either commonly 
or differentially. This curvature radius error corresponds 
to the following mirror figure error: 



Az(x,y) = 1.04nm 



-1 



6cm 



y 

6cm 



(26) 

where Az(x,y) is the mirror surface height error, and 
Az(x,0) is show in Fig.0 

We show the loss of the SNR in Fig. [3J With changes 
of the SRC degeneracy, characterized by the Guoy phase 
and the SRC g-factor, we see significant change in the 
loss of the SNR. In the baseline degenerate design, for 
common and differential perturbations we lose 4% and 
0.4% of the SNR. Note that at leading order the SNR loss 
is proportional to the square of the size of the error, so 
we have, for instance, 4 times the above SNR losses when 
Ai?iTM = 10m instead of 5m. The bigger of these SNR 
losses is consistent with the estimate based on geometric 
optics approximation in Section IV J of Ref. |l3| . 

Figure shows that, when the degeneracy is reduced, 
the SNR loss drops by more than two orders of magni- 
tude, making the curvature error harmless. The most 
striking features in the plots are the peaks corresponding 
to huge SNR loss at some nondegenerate SRC configura- 
tion. This happens when the Guoy phase of the HOMs 
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FIG. 2: Mirror figure errors given in Eqs. 1261 and 1271 . We 
show as solid curves the error Az in units of nanometers along 
the a;-axis [i.e., Az(a;,0)]. The dashed curve is the power 
profile of the fundamental Gaussian mode plotted in arbitrary 
units. The beam size at the mirror is w — 6cm (c.f. Eq. 10), 
so the variance of the power profile is a — 3cm. The peak- 
to-valley figure errors inside the 2er (6cm) region are about 
2.0nm and 9.3nm for the two cases. 



eter, and when the SRC is degenerate, because the carrier 
light is not anti-resonant in it, it behaves like "resonant 
carrier extraction," and sucks carrier light power out of 
the AC thus reducing its power build-up. When the SRC 
is nondegenerate, HOMs are rejected by both the SRC 
and the AC, and the carrier power builds up as usual in 
the AC in almost entirely its fundamental mode, without 
losing any significant power. The SNR loss in this case 
is due to the loss of carrier power in the AC, since the 
signal sideband HOMs are coupled to the symmetric port 
and rejected by the PRC, which is not properly detuned 
for the signal sideband frequency. 

In conclusion, common and differential errors in the 
ACs reduce the SNR in different ways, and the dif- 
ferential errors reduce the carrier build-up significantly. 
Finally, as indicated above by "resonant carrier ex- 
traction," the differential curvature error sends a huge 
amount of carrier power in HOMs toward the dark port 
(7.5W in HOMs and 0.3W in fundamental modes inside 
the SRC for Ai?iTM = 5m), and thus sends reference 
light that is mostly in HOMs toward the photon detec- 
tor. This HOM light must be cleaned out by an output 
mode cleaner. 



(in this case, (2,0) and (0,2) Hermite-Gaussian modes) 
cancels the SRC detuning (<fif IC = 0.06 for broadband 
design), so that the HOMs of both the carrier and the 
signal light are resonant in the SRC while the fundamen- 
tal modes are detuned. This is clearly a bad choice of 
SRC degeneracy. We refer to it as the HOM resonant 
peak. When more HOMs are coupled into the interfer- 
ometer by perturbations, we should avoid all such cavity 
configurations in which some HOM (of an order not so 
high that it suffers strong diffraction loss) has a Guoy 
phase ?7hom that nearly cancels the SRC tuning phase. 

Figure0]shows the effect of the SRC degeneracy on the 
loss of carrier light due to differential curvature errors 
on the ITMs. In addition to the HOM resonant peaks 
discussed above, there are other noticeable features in 
the carrier light. In the AC, curvature errors reduce the 
carrier-power build-up in ways described below. When 
the SRC becomes nondegenerate, for common and dif- 
ferential curvature errors on the ITMs, the carrier and 
signal sideband light behaves very differently (compare 
Figs. 0and EI). 

For common curvature errors, the HOMs are coupled 
into the symmetric port of the interferometer, which is 
accepted by the degenerate PRC but are anti-resonant in 
it, so they get reflected back into the AC. Since the HOMs 
do not enter the SRC, its degeneracy causes no difference. 
In this case the carrier light power is controlled by the 
AC and is hardly affected by the mirror deformations, 
thus we do not show it. The SNR loss in this case is due 
directly to the coupling of signal-light HOMs to the dark 
port. 

For differential curvature error, by contrast, the HOMs 
are coupled to the anti-symmetric port of the intcrfcrom- 



B. Different modes of deformation on the SRM 
and ITMs 

In this section, we consider ACs and the SRC with 
different modes of deformation in the SRM and ITMs. 
In Fig. [3 we show the loss of the SNR as a function of 
SRC degeneracy level, when there are common curvature 
radius errors of AR = lm on the ITMs, and a fourth- 
order polynomial deformation on the SRM with the form 



Az(x,y) 



-4.65nm+ 17.94nm 

L ( — 

\6cm 



.6cm 



6cm 



V 



6cm 



(27) 



This Az is plotted in Fig. HJ 

Comparing Fig. [S] with Fig. El we see two more HOM 
resonant peaks generated by the fourth-order Hermite- 
Gaussian modes at Guoy phase around 0.8 and 1.4 ra- 
dian. We see from Fig. [5] that, although the optical 
modes in the AC and the SRC are different, and there 
is an eigenmode mismatch on the ITMs, the nondegen- 
erate SRC does not simply reject the mismatched part 
of the signal power, and the loss of the SNR is still 
very low (away from the HOM resonant peaks). This 
is because the SRC effectively reflects the HOMs that 
cause the mode mismatch back into the AC and helps 
the fundamental mode build-up in the AC, so long as 
the SRC tuning for the fundamental mode resonance is 
unchanged. (This tuning depends on how the control sys- 
tem works. Modeling of the control sidebands is being 
considered in more sophisticated simulations that are un- 
der development, e.g., Advanced-LIGO FFT simulations 
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FIG. 3: Loss of the SNR in Advanced-LIGO interferometers due to mirror curvature radius error on the ITMs, at the resonant 
signal sideband frequency fo — X/2n, as a function of the SRC degeneracy level. The curvature radius of the ITMs is Ritm = 
2076.4m and we consider an error ARitm = 5m, which is equivalent to the mirror surface height error given in Eq. (1 "_2 ( > t . We 
consider the two cases where the curvature errors on the two ITMs are common (i.e., curvature radii are -Ritm — ARitm) and 
differential (i.e., curvature radii are -Ritm ± ARjtm), and show the loss of the SNR in the left and right panels, respectively. 
The horizontal axis is the level of degeneracy, measured by two quantities: the g-factor of the SRC and the one-way Guoy 
phase in the SRC. The horizontal position of the left-most point of the curve corresponds to the degeneracy level of the SRC 
in the current Advanced-LIGO baseline design. 




FIG. 4: Change of the carrier light power in the AC and SRC in Advanced-LIGO interferometers due to differential mirror 
curvature radius error on the ITMs, as a function of SRC degeneracy. We consider the same differential curvature error as in 
Fig. [3] In the left panel, we show the loss of carrier power (all optical modes) in the AC in units of kW (the ideal AC carrier 
power is about 825k W); in the right panel we show the increase of the carrier power (all optical modes) in the SRC in units 
of W (the ideal SRC carrier power is about 1W). The horizontal axis is the level of degeneracy, measured by two quantities: 
the g-factor of the SRC and the one-way Guoy phase in the SRC. The horizontal position of the left-most point of the curve 
corresponds to the degeneracy level of the SRC in the baseline Advanced-LIGO design. 



|21|). More precisely, this argument is valid because the 
coupling of the carrier to the signal sideband happens 
at the ETM, outside the SRC; if it happened inside the 
SRC, the non-matching HOMs would be expelled and 
thus the signal power in the fundamental mode would be 
reduced. 



C. Curvature radius error on the ITMs: 
Narrowband configuration 



All examples above are simulations for the broadband 
Advanced-LIGO configuration, i.e., the RSE configura- 
tion, in which the signal storage time in the AC is re- 
duced by the SRC. In the narrowband configuration, by 
contrast, the signal light is truly being "recycled," and 
the storage time in the SRC is an order of magnitude 
longer than in the RSE scheme, which could change our 
broadband-interferometer results significantly. In Fig. 
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we show simulation results for the narrowband configu- 
ration with differential curvature errors AR = 2m on the 
ITMs. The SNR loss for the baseline degenerate SRC 
is about 5%, which is very significant. In Ref. [l^, the 
narrowband configuration also implies the most severe 
constraint on mirror figure errors. The SRC finesse in 
the narrowband configuration is much higher than in the 
broadband Advanced-LIGO configuration, so the HOMs 
being excited inside the SRC are built up to higher power, 
when the SRC is degenerate. As for other configurations, 
the loss of the SNR is reduced by orders of magnitude 
when the SRC is changed from degenerate to nondegen- 
erate. 



FIG. 5: Loss of the SNR in Advanced-LIGO interferometers 
due to mirror curvature radius error on the ITMs and the 
SRM, and some higher-order mode deformation on the SRM, 
at the resonant signal sideband frequency /o — A/2-7T, as a func- 
tion of SRC degeneracy. We consider a common curvature 
radius error on the ITMs with A7?itm = lm, i.e., one- fifth 
the error considered in Fig. and 0] or equivalently one- fifth 
the mirror surface height error given by Eq. fRft. The mirror 
surface height error of the SRM is given in Eq. I|27|l . The 
horizontal axis is the level of degeneracy, measured by two 
quantities: the g-factor of the SRC and the one-way Guoy 
phase in the SRC. The horizontal position of the left-most 
point of the curve corresponds to the degeneracy level of the 
SRC in the baseline Advanced-LIGO design. 
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FIG. 6: Loss of the SNR in narrowband (see Eqs. If 8H and 
1191 for parameters) Advanced-LIGO interferometers due to 
common mirror curvature radius errors on the ITMs, at the 
resonant signal sideband frequency /o — X/2n, as a function 
of SRC degeneracy. The common error of curvature radius is 
A_Ritm = 2m, i.e., two-fifth the error considered in Figs. [3] 
and 2] or equivalently two- fifth the mirror surface height error 
given m Eq. The horizontal axis is the level of degen- 

eracy, measured by two quantities: the g-factor of the SRC 
and the one-way Guoy phase in the SRC. The horizontal po- 
sition of the left-most point of the curve corresponds to the 
degeneracy level of the SRC in the current Advanced-LIGO 
baseline design. 



D. Optimal SRC degeneracy 

From the examples above, we see that the mode mix- 
ing and consequent problems are suppressed by making 
the SRC nondegenerate. The optimal Guoy phase in the 
SRC for Hermite-Gaussian modes is the range 0.2 to 1.3 
radians. In the examples above we showed the loss of 
the SNR only at the most sensitive signal sideband fre- 
quency [i.e., f - X g /2n given by Eq. P0|l]. In Fig. [7| 
we show the interferometer response of various signal fre- 
quencies assuming that the SRC is either degenerate with 
the current design parameters or nondegenerate with the 
MMT design suggested by Muller 01 ( tne cavity Guoy 
phase is 0.38 radians, and thus nondegenerate). There 
is a strong suppression of the SNR loss across the en- 
tire Advanced-LIGO sensitive band, and a small shift of 
the most sensitive frequency (remember however that at 
low frequency / < fOOHz, we should also consider the 
radiation pressure noise) . 

At first sight, there seems to be a wide optimal range 
from which to choose the SRC degeneracy, but our free- 
dom is actually quite limited. One obvious constraint is 
that we need to avoid those values of the Guoy phase 
that give rise to the HOM resonant peaks. In a realis- 
tic interferometer, with many more HOM perturbations, 
there would be a large number of HOM resonant peaks 
across the optimal Guoy phase range, so the degeneracy 
should be chosen carefully, through careful simulations. 

However, the worst difficulty posed by the above Guoy 
phase range is a practical one. There are two obvious 
ways to achieve the desired g-factor: reduce the beam 
size in the SRC, or increase the length of the SRC. Un- 
fortunately, to achieve the required low degeneracy, we 
need cither a very small beam waist size near the SRM, 
or a kilometers-long recycling cavity length. In the next 
section, we discuss the practical alternative designs of 
these two types. 

V. ALTERNATIVE DESIGNS 

There are two obvious ways to reduce the SRC degen- 
eracy: reduce the beam size and/or increase the optical 
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FIG. 7: Loss of the SNR in an Advanced-LIGO interferometer due to mirror curvature radius error on the ITMs, for signal 
sidebands with various frequencies. The mirror curvature errors are the same as those consider in Fig. [3] and we show in the 
left and right panels the results for common and differential errors. In each case, we consider two designs of the SRC: one 
degenerate SRC with the Advanced-LIGO baseline design parameters, and one nondegenerate SRC with one-way Guoy phase 
77 = 0.38 (i.e., g-factor 0.86) prescribed preliminarily in the MMT design [TUl . In each plot the solid curve with larger variations 
shows the loss of the SNR when the degenerate SRC is used, while the solid curve with smaller variations shows the loss of the 
SNR when the nondegenerate SRC is used. The dashed curve, plotted in arbitrary units, shows the amplitude of the output 
signal sideband light in an ideal Advanced-LIGO interferometer with fixed GW strain. The maximum of the dashed curve is 
at the frequency fo — \ g /2n with X g given by Eq. 1201 . 



path length. We discuss two ideas in turn. 

To reduce the SRC beam size, we could add a lens in 
the recycling cavity; but to get a Guoy phase between 0.2 
and 1.3 radians, the beam must be focused so strongly 
that the waist size of the beam is of order 3 x 10 _5 m, 
and the waist location must be tuned precisely to a few 
millimeters away from the SRM. For a 1 x 10 _5 m beam 
size on the SRM, we will have a 10GW/m 2 power density 
heating one spot on the mirror, which is too high to be 
practical. This problem was pointed out qualitatively by 
Bochner in his FFT simulation work. 

We can circumvent this problem by introducing multi- 
ple steps of beam focusing, i.e., bringing down the beam 
size step by step with more optical elements, so there 
is substantial Guoy phase accumulated during each step 
with small beam size. It is clearly preferable to use reflec- 
tive optical elements in this scheme. A practical design 
based on moving the mode-matching telescopes (MMT) 
into the recycling cavities has been proposed by Muller 
and Wise in Ref. [ll| . More specifically, the MMTs out- 
side the recycling cavities used to match the beam size 
from mm scale to cm scale between the light source and 
the PRM, and from cm scale to mm scale between the 
SRM and the photo detector are moved into the cav- 
ities to reduce the beam size there. Two mirrors are 
used in each MMT. The first brings down the beam size 
from ~ 6cm to millimeter scale, and the second tunes 
the shape of the millimeter-scale beam to achieve the de- 
sired degeneracy. Practical parameters for the PRC can 
be found in Ref. and the coupled recycling cavities 
and cavities in the MMT arc stable in principle. The 
MMT introduces more mirrors and cavities into the in- 
terferometer, and there are more control problems for 



the new mirrors. Experimental studies of issues about 
implementing this MMT design are ongoing. 

Another way to reduce the SRC degeneracy is to use a 
longer SRC, though the length must be on the kilometer 
scale to achieve the necessary g-factor. A natural idea 
is to bend the SRC into one of the beam tubes of the 
ACs, and make it 4km long. This design seems to have 
the immediate problem of light scattering noise in the 
crowded arm tubes, but it does have important scientific 
advantages. The 4km SRC idea has been advocated for 
a long time, and reducing the SRC degeneracy is just 
one gain among other advantages. Mizuno |(| suggested 
using a 4km long SRC to collect power in both signal 
sidebands and increase the SNR by a factor of 2. Buo- 
nanno and Chen ,8], in their analysis of beating the SQL 
with a signal-recycled interferometer, found that the gain 
in peak sensitivity is vulnerable to optical losses in the 
short SRC, and km long SRC with fewer light bounces 
might solve this problem. Moreover, a long SRC intro- 
duces a frequency-dependent correlation between the two 
quadratures of the vacuum field, and might thereby bring 
interesting changes to the optical noise spectrum of the 
interferometer. Finally, light scattering in the beam tube 
may not cause any significant problems according to pre- 
liminary estimates, and more investigation of this issue 
is ongoing. 



VI. CONCLUSIONS 

We have set up a simplified model of an Advanced- 
LIGO interferometer, with mirror deformations in the 
form of figure errors, and we have simulated the carrier 
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and signal sideband optical fields in this model interfer- 
ometer by a mode decomposition method. Using our 
simulations, we have investigated the loss of the SNR 
when SRCs with various degeneracies are used, and we 
have found an optimal range of SRC degeneracies that 
minimizes the loss of the SNR. 

For the current degenerate SRC design, we found that 
the loss of the SNR, due to mode mixing between the 
fundamental mode and HOMs, is severe, and our results 
are consistent with the geometric optics estimates made 
in Ref. 01 . Assuming that the loss of the SNR scales 
quadratically with the size of the mirror figure error, and 
requiring the SNR loss due to mirror deformation to be 
smaller than 1%, we found a very severe constraint on 
the mirror deformations. In the broadband Advanced- 
LIGO design, If the only type of mirror deformation is 
curvature-radius error, this error must be smaller than 
2.5m-7m on the ITMs, which corresponds to peak-to- 
vallcy figure errors smaller than lnm-3nm in the central 
region of the mirrors. In reality, the mirror deformation 
is formed by a combination of many spatial modes, and 
when we consider the next lowest order spatial modes, as 
depicted in Eg. 1271 and Fig. [3 we get a constraint ~ 4nm 
on the SRM. In the narrowband design, the constraint on 
the ITMs is tightened to ~ 0.4nm. Considering the fact 
that at leading order, the loss of the SNR grows quadrat- 
ically with the size of the figure errors, and losses due to 
figure errors with different spatial modes are added lin- 
early, the constraint on mirror figure errors and thermal 
effects is very difficult to achieve with current technol- 
ogy. Another minor problem we found with the baseline 
degenerate SRC is that, when there are differential per- 
turbations on the ACs, a huge amount of carrier light in 
HOMs is coupled to the dark output port and overwhelms 
the fundamental- mode reference light. To remove this 
large HOM light at the output requires a reliable output 
mode cleaner. 

We have shown that a nondegenerate SRC could solve 
this problem, by suppressing the mode mixing and re- 
ducing the loss of the SNR by orders of magnitude, in 
the Advanced-LIGO sensitive band. We have shown also 
that a nondegenerate SRC does not simply reject mode 
mismatched light from the AC; it also helps the funda- 
mental mode to build up. 

We propose using a nondegenerate SRC in Advanced 
LIGO with the one-way Guoy phase between 0.2 and 1.3 
radians. There are difficulties in achieving this optimal 
degeneracy in practice, and we have discussed two pos- 
sible alternative designs for doing so. In the first design, 
we move the MMT into the recycling cavities and in the 
second, we use a 4km long SRC. 

A more complete simulation of the optical fields in- 
side an Advanced-LIGO interferometer u sing FFT prop- 
agation methods is under development |2l|. This will 
effectively include hundreds of HOMs, and will model 
important physical factors such as thermal effects on the 
mirrors, and control sideband fields. This new simulation 
is aimed at mapping out the phase fronts of the light in 



a very realistic Advanced-LIGO model, to an accuracy 
of ~ 10~ 6 . This new simulation, among other goals, will 
help perfect designs of SRCs with optimal degeneracy. 
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APPENDIX A: ABBREVIATIONS AND 
SYMBOLS 

There is a large number of abbreviations and symbols 
used in this chapter. We give a full list of them in this 
appendix, for easy reference. 

Abbreviations: 

AC arm cavity 

BS beam-splitter 

EM electromagnetic 

ETM end test mass 

FFT fast Fourier transform 

GW gravitational wave 

HOM higher-order mode 

ITM input test mass 

MMT mode-matching telescope 

PRC power-recycling cavity 

PRM power-recycling mirror 

RF radio frequency 

RSE resonant sideband extraction 

SNR signal-to-noise ratio 

SRC signal-recycling cavity 

SRM signal-recycling mirror 

TCS thermal compensation system 
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Symbols: 

Az mirror surface height error 

r\ Guoy phase 

A carrier light wavelength 

X g resonant signal sideband frequency in the 

combined SRC- AC cavity 
ljq carrier light frequency 
uj g GW sideband frequency 
<j) a phase shift in cavity "a" 
4>h phase shift due to GW strain 
c speed of light in vacuum 
d Michelson arm length difference 
/o carrier light frequency 
k carrier light wave number 
L Length of arm cavity 
l a common length of recycling cavity "a" 
M a reflection operator of mirror "a" or 

of cavity "alpha" as a compound mirror; 

Exceptions: Michelson operators Mq , M§, 
and Mg defined in Eq. (|U9|) 
P a one-way propagator in cavity "a" or 

propagator of path "alpha" 
r a amplitude transmissivity of mirror "a" 
T a transmission operator of mirror "a" 
t a amplitude transmissivity of mirror "a" 
Z distortion function 



Superscripts: 

B broadband design 
C carrier light 
N narrowband design 
S signal sideband 



APPENDIX B: ASSUMPTIONS AND 
APPROXIMATIONS IN ADVANCED-LIGO 
INTERFEROMETER MODEL 

In this appendix, We list the principal assumptions 
and approximations adopted in our simplified Advanced- 
LIGO interferometer model, that are introduced and dis- 
cussed in various places in the chapter: 

(1) We consider only the five lowest-order HOMs and 
focus mostly on the two lowest-order modes, i.e., modes 
excited by mirror curvature-radius errors 

(2) We assume the sizes of mirrors in the interferome- 
ters are large enough for diffraction losses to be negligible. 

(3) We assume that the interferometer noise is dom- 
inated by the photon shot noise, and ignore radiation- 
pressure noise. 

(4) We assume all degrees of freedom of the interfer- 
ometer (cavity tuning, optical clement alignment, etc.) 
are fixed to their ideal values as if there were no mirror 
deformation, except for the tuning of the PRC and the 
AC. The PRC and AC are tuned to maximize the total 
carrier light power (sum of power in all optical modes) 
in the AC, for fixed input light power. 

(5) We include in our simulation only the carrier light 
and one signal sideband, without any sidebands for con- 
trol purposes. 

(6) We choose a microscopic asymmetry between the 
Michelson arm lengths that sends about 1W of carrier 
power to the SRM. 

(7) We assume that the beam-splitter is perfect. 

(8) We adopt the short distance approximation 
(Sec. IHH. 

(9) We ignore the difference in phase shifts in the short 
recycling cavities between the carrier and signal side- 
bands. 

(10) We assume that all mirrors are lossless. 



Subscripts: 

acl online arm cavity 
ac2 offline arm cavity 
b beam-splitter 

C common mode in Michelson arms 
D differential mode in Michelson arms 
±d a ±d trip in one of the Michelson arms 
el end test mass of online arm cavity 
c2 end test mass of offline arm cavity 

11 input test mass of online arm cavity 

12 input test mass of offline arm cavity 
p power-recycling mirror 

pre power-recycling cavity 

rt acl round trip in online arm cavity 

rt ac2 round trip in offline arm cavity 

s signal-recycling mirror 

sre signal-recycling cavity 



APPENDIX C: SOLVING FOR THE OPTICAL 
FIELDS IN THE INTERFEROMETER 

In this appendix, we write out the explicit form of the 
interaction operator M in the Hermite-Gaussian modal 
space, and write equations in the form of Eq. 112 111 for all 
the carrier and signal sideband fields displayed in Fig. ^ 
We solve these equations analytically in terms of the in- 
put field £?inOj the GW strain, and the propagation and 
interaction operators. The SNR is proportional to the 
fundamental mode component of the output signal side- 
band field E^ ut 

The carrier light and the signal light have inputs at 
different positions (the signal input is effectively at the 
ETMs), so their solutions are different and have to be 
treated separately. In the following we will denote with 
superscripts "C" for carrier and "S" for signal. Further- 
more, because of their different frequencies, their propa- 
gation operators have to be distinguished as well in the 
AC; however, for the short recycling cavities, since the 
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fractional difference between the carrier and signal fre- 
quencies is on the order of 10~ n , the phase difference 
between carrier and signal light is negligible and we use 
the same SRC operators. For the same reason, the inter- 
action operators are the same for carrier and signal light, 
since the perturbation effect has effectively a length scale 
~ lOnm. 

For Hermitc-Gaussian modes, the matrix element of 
A I in the modal space is given in Ref. ^3 as 



mn\e 
mn 



expf-i/c ^2 \°P) z o P ,qr(qr\ 



kl 



(CI) 

where Z op ^ qr can be calculated through the following 
quantities: 

(i) Cij is given by the expansion of the distortion func- 
tion Z{x,y) into Hermite polynomial modes: 



-kZ(x,y) = ^2 c iJ H i 

i,j V 



2x\ T ( V2x 



w{z) 



w{z) 



(C2) 



in which w(z) is the beam size at the reference surface 
where the operator is defined; 

(ii) h]? t qr is given by the expansion of Hermite polyno- 
mial products: 

Hi (x)Hj (y)H q (x)H r (y) = £ h% qr H s (x)H t (y) . (C3) 



(iii) T l J v qr is given by the integral 



7*j 

op,qr 



{op\H. l {x)H J {y)\qr)\ z=() 



(C4) 



which is given analytically when the mirror size is infinite 
as 



rplj 

op,qr 



,L st,qr 



' 2°o\2Pp\ 
2iq\2 r r\ 



(C5) 



When the mirror size is finite, T l J v qr has to be calculated 



numerically. 

(iv) Z op ^ qr is then given by: 



J op,qr 



"cut* 



i,3 



(C6) 



Therefore, given Z(x,y), which embodies information 
about the mirror distortion, and the beam size at the 
optical clement w(z), we get the representation of the op- 
erators M mrit ki in modal space following the steps above. 

Now we turn to the optical fields. As described in 
Sec. IIIII we solve for the carrier field as the first step. 
Using results for LIGO that have been derived in Ref. [l7j 
and Ref. (2^, and following the conventions in Sec. IIIII 
we define the following. 



(i) The round-trip propagator in the AC is given by 



frtacl = mrelMiiP^MerP^ 
^r?ac2 - r i2 r c2 M {2 Pg 2 M e2 Pg 2 



(C7) 



again we use subscripts "1" and "2" to denote the online 
and offline arm cavities. 

(ii) The reflection operators of the ACs (for reflecting 
off the AC from the ITM side) are given by 

A4 C c 2 = n 2 (m^-^TuM^P^I-P^)- 1 ^ , 

(C8) 

where M[ x and M[ 2 are reflection operators of the ITMs 
as seen from the recycling cavity side. 

(iii) The Michclson cavity operators are given by 

MS = tlP<tM£ 1 P<t + rlP° d Mg 2 PS d , 

Mg = t hs r hs (P°MC cl pC + PE d Mg 2 PS d ) , (C9) 

where P d and P G d are propagators through the differen- 
tial length d of the Michelson arms [Eq. Ijllljl ] . By com- 
bined with Pp IC and P G C , i.e., the propagators through 
the common lengths of the PRC (Z prc = 8.34m) and 
SRC (Z src = 8.327m), we get the propagators from the 
PRM to the PRM and the SRM (i.e., the "common" 
and "differential" mode propagators): P p c [c MpF p c rc and 
pC m g P g 

The coupled equations for the carrier fields are then 

pC _ j- pC pC pC _ pC pC pC 
■C's — t bs^p rc ^d "rel ' bs^p lc ^_d- C/ re2 J 

E c = r\, s P G c P d E G cl + t]j S P G c P G d E G 2 , 
E° = -r p M p E G + t p T p E in0 , 
Eg = -r s M s £ a c , 

^inl = *bs-Pprc-fd ^si + r bsP£ c Pd ^ar ) 

pC _ _ pC pC pC i j- pC pC pC 
Mn2 — ' bs^prc-T-d^sr ^ ( 'bs-r src -r_d J ^ar ' 

E? e i = M G cl E G nl , 



Ere2 — M ac2 E- m2 ) 



(CIO) 



and the circulating fields inside the AC are given in terms 
of Eg, and by 



^cirl — ^il-Pacl [J •f'rtacl) ^il-^inl i 

Eg 2 = k 2 Pg a (I - P^Y 1 T i2 Eg 2 . (Cll) 
Solving Eq. 1C10|) , we have 

E° = tJl + r p M p P G M%P° - r p r s M p P G c M§Pg 



30 



x(I + r s M s p£ c M§P° IC ) -M S P£ C M£P£ C ) T p E in0 , 



C A/fC nC 



sre 

1 



(J + r s M s pC c M§p£ c ) 1 r s M s P s c IC M£p G c E G 



(C12) 



15 



All other carrier fields can be easily calculated from E^. 
and E°. 

Assuming a monochromatic GW wave passing through 
the interferometer and shaking the ETMs differentially 
with strain ho cos u) g t, the carrier light that is incident on 
the ETMs (i.e., £*? rl and E^ ir2 ) is coupled to the motion 
of the ETMs and generates signal sidebands at frequen- 
cies u> ± ojg with the EM fields given by 

£ s S i g i = i<t>hr e iM e iE° hl , 

£ s S ig2 = ^ h r c2 M c2 Eg. 2 , (C13) 

where 4>h = 2kh L is the phase shift due to the GW 
strain. These fields are the input for the signal-light field 
in the interferometer. 

For the signal sideband field, we define the round-trip 
propagator of the AC: 

^ S tad = n 1 r el M el Pg x M il Pg l , 

^ S tac 2 = n 2 r c2 M c2 Pg 2 M- l2 Pg 2 , (C14) 

and we define for the signal sideband the reflection op- 
erators of the ACs Mf cl and M^ c2 , and Michelson cavity 
operators: A'Iq and Mq, in the same way as their carrier- 
light counterparts were defined in Eq. I|C9J) . but using 
signal sideband propagators. 

The coupled equations for the signal fields arc similar 
to Eq. I|C10(I : we only need to change all quantities for 
the carrier fields in Eq. I|C10|) to their counterparts for 
the signal fields, and change the positions of the input 
fields in three of the equations: 

£ s s r = -r p M p £ s s , 

^rel = ^acl-^inl + ^U^il-Pacl (l ~ ^acl) ^sigl ! 
E ro 2 = M^. 2 E? n2 + t- l2 T l2 P^ c2 (/ - Pfta.c 2 ) E aig 2 ■ 

(C15) 



Solving Eq. (|Ul5)l for we obtain the output signal 
field: 



P D s ut = t s T s E s a = t s T s (I + r s pS c M*pS c M s ) 

x (-r^M^M^ + P s s iga ) , (C16) 



where 

Ef = (/ + r p P p s rc Af s p s cMp _ rpTaP S c A/SpS cMs 

x (/ + rsP^MSp^Msy 1 P^MgpS^/p) ~* 
x (ps gs _ r a P^M^ IC M a 

x (I + r 8 P 8 s rc MgP s s rc M 8 )- 1 P s s iga ) , 

-^edgs = ^il^bsPprc-Pd ^il^acl (j — ^it acl) ^sigl 

— ^2f*bsPp rc -P-d^ 1 i2-Pac2 (l — ^rtac2) ^sig2 > 

-^figa = iilTbs-Psrc-Pd ^il-Facl " ^rtacl) ^figl 

+ti2ib S P S re-P-d^i2-Pac2 — ^rtac2) -^fig2 ■ 

(C17) 

According to Eq. (J2BJ), the SNR is proportional to the 
amplitude of the output signal sideband field E^ nt in 
the fundamental optical mode. The fractional loss of 
the SNR is then the fractional loss of the fundamental 
mode amplitude in E^ ut when mirror deformations are 
introduced through the mirror reflection and transmis- 
sion operators M a and T a . 
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